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Abstract 

Scanning tunneling microscopy and spectroscopy in magnetic field was used to study Landau 
quantization in graphene and its dependence on charge carrier density. Measurements were carried 
out on exfoliated graphene samples deposited on a chlorinated SiC>2 thermal oxide which allowed 
observing the Landau level sequences characteristic of single layer graphene while tuning the density 
through the Si backgate. Upon changing the carrier density we find abrupt jumps in the Fermi 
level after each Landau level is filled. Moreover, the Landau level spacing shows a marked increase 
at low doping levels, consistent with an interaction-induced renormalization of the Dirac cone. 



PACS numbers: 73.22.Pr, 71.70.Di, 73.20.At, 73.43.Jn 



One of the hallmarks of the relativistic charge carriers^ 2 in graphene is the appearance 
in a magnetic field of an unusual Landau level (LL) at zero energy which reflects the chiral 
symmetry of the low lying excitations. The presence of this LL has been inferred in magneto- 
transport measurements employing the standard configuration of graphene supported on 
SiC>2 3 4 from the conspicuous absence of a quantum Hall plateau at zero filling-factor. Re- 
markably because in graphene the carriers reside right at the surface, the LLs (including the 
LL at zero-energy) can be accessed directly through scanning tunneling spectroscopy (STS) 
as was demonstrated in studies of graphene samples supported on the surface of graphite^. 
However, the LLs were not observed in STS measurements on graphene samples supported 
on insulating substrates which allow control of the carrier density through gating. This is 
because due to the purely two dimensional nature of graphene, substrate induced potential 
fluctuations obscure the intrinsic physics of the charge carriers close to the Dirac point. 
One way to overcome this limitation is to use suspended samples^ where transport mea- 
surements have shown that in the absence of the substrate the intrinsic Dirac point physics 
including interaction effects is revealecP^. The use of suspended samples is however limited 
due to their fragility, small size and reduced range of gating. Finding a minimally invasive 
insulating substrate on which graphene can be gated and also visualized is therefore of great 
interest. 

By using scanning tunneling microscopy (STM) and spectroscopy (STS) we show that 
for graphene supported on Si02 substrates which were treated by chlorination to minimize 
trapped charges and in sufficiently large magnetic fields, the LL sequence specific to single 
layer graphene and its dependence on carrier density can be accessed. Upon varying the 
carrier-density sudden jumps of the Fermi-energy are observed after filling each LL. Moreover 
the measured density-dependence of the LL spacing shows a rapid increase upon approaching 
the Dirac-point, consistent with an interaction-induced renormalization of the Dirac cone. 

A simple method for preparing high quality graphene samples is mechanical exfoliation 
from graphite followed by deposition on the surface of Si02 capping a Si crystaP. This 
relative ease of sample preparation has stimulated numerous studies aimed at probing the 
extraordinary properties of graphene which range from novel physical properties to a wide 
range of applications. A similar method was used for the graphene devices studied here. 
Graphene samples prepared by exfoliation from HOPG graphite were deposited on the sur- 
face of a p-doped Si wafer capped with 285 nm of Si02- To minimize substrate induced 
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FIG. 1: (a) Sketch of the experimental set up showing a graphene flake deposited on S1/S1O2, 
the STM tip above and the back gate connected to the graphene. (b) STM topography image of 
a typical 300 nm x 300 nm graphene area. Tunneling conditions: It = 20 pA, Vbi as = 190 mV.(c) 
STM atomic resolution image on the graphene flake (It = 20 pA, Vuas = 300 mV).(d) Profile 
line through the topography image (at the position indicated in (b)) showing the corrugation of 
graphene on Si02 substrate. 

disorder we used a chlorinated thermal oxide^ purchased from Nova Electronics which was 
annealed in forming gas at 230 °C prior to the deposition of graphene. The metallic leads 
were patterned using e-beam lithography followed by evaporation of 1.5/30-40 nm Ti/Au. 
Prior to measurement, the samples were baked overnight in forming gas at 230 °C. STM and 
STS measurements were performed at 4.4 K in a home built STM (Fig. [T|(a)). After baking, 
the surface of graphene is found to be free of fabrication residues on areas as large as 1 /zm 
x 1 fim as measured by low temperature STM. A typical 300 nm x 300 nm area is shown 
in Fig. [T|(b) . As found previously the samples are rippled on a fine scale of 2-5 A in the 



vertical direction and a few nm laterally 22113 . When zooming into atomic resolution (Fig. 
[ljc)) the typical honeycomb structure observed is a first indication of the sample cleanliness 
and quality. 

The intrinsic density of states (DOS) of graphene in the absence of field is known to be 
V-shaped and to vanish at the Dirac point. In samples where substrate induced disorder 
is negligibly small: graphene flakes on graphite, it was possible to measure this intrinsic 
spectrum in zero fielcP. However, when graphene is deposited on an insulating substrate 
such as Si0 2 its properties are modified due to defects, strain, inhomogeneous doping and 
trapped charges 23 ^. These altering effects are responsible for a wide range of non-intrinsic 
spectra reported in the literature for graphene samples deposited on SiC>2 substrate j 16 ' 17 1 19 1 . 
In our case, scanning tunneling spectra in the absence of field and for a neutral gate show 
two dips, one at the Fermi energy (Ep) where E=0 eV, and the other at w 170 meV (Fig. 
(2^c)). The deviation from the expected V-shape is attributed to substrate-induced disorder 
which smears the Dirac point. Similar to results reported m 1 15 * 19 ! but unlike those irP^, we 
do not see evidence for a the opening of a gap at E^. 

The inhomogeneous carrier density seen as electron-hole puddles close to the Dirac Point 
(DP)IMI£I is clearly visible in the local density of state map shown in Fig. [2^b) that corre- 
sponds to the region whose topography is shown in Fig. [2^a). We estimate a typical puddle 
size of d=20 nm. This imposes a length scale and a corresponding energy scale which 
separates between disorder-controlled and intrinsic phenomena. In finite field the intrinsic 
physics of the charge carriers will thus become apparent only when the cyclotron orbit is 
smaller than the characteristic puddle size, l c = = 25.64 nm/y/B m d/2, where B is 
the magnetic field normal to the layer. This defines a characteristic field, B c = below 
which the LL spectrum is smeared out by disorder. Consistent with this picture we find 
no distinct LL features in the tunneling spectrum for fields below B c 6 T (Fig. [3[a)), 
corresponding to l c (6 T)=10.5 nm w d/2. Above 6 T the spectra develop peaks that evolve 
and become more pronounced with increasing field as shown in Fig. (3|b). To compare the 
field and level index dependence of the measured spectra to that expected for the quantized 
LL energies specific to the two dimensional Dirac electrons in graphene: 

E = ±v F ^2eh \N\ ■ B, N = 0, ±1, ±2, ... (1) 

we plot in Fig. [3|c) the measured peak energies as a function of the reduced parameter 
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FIG. 2: (a) Topography image corresponding to the Local Density Of States (LDOS) map in (b) 
(It = 20 pA, Vbias = 140 mV). (b) LDOS map taken at E=140 meV showing doping inhomogeneity 
in graphene. (c) STS spectrum at B=0T showing smearing of the Dirac Point (indicated by ED) 
due to disorder. 

sgn(N) ■ \/\N\ B. We find a linear dependence in agreement with the expected LL sequence 
for single layer graphene and from the linear fit we extract the Fermi velocity: vp=(1.07 
± 0.02) x 10 6 m/s. We find that there is a 5-10% variation in the Fermi velocity across 
the sample. This value is consistent with other measurements of the Fermi velocity in 
graphene on SiO j 16 * 20 ^. Compared to the results obtained for graphene on the surface of 
graphite^ and for epitaxial graphene on SiCP^, the levels are broader and their width varies 
with position suggesting a shorter carrier lifetime due to extrinsic scattering mechanisms 
such as trapped charges, ripples, etc. The typical line-width for LLs with N^0, 20-30 meV, 
corresponds to carrier lifetimes, r « xjfjT ~ 22 fs - 32 fs consistent with values found by 
other methods in non-suspended graphene samples^ 18 * 20 * 21 * 23 ! 
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FIG. 3: Landau levels in ungated graphene on SiO"2 (Vg=0 V). (a) Scanning Tunneling Spec- 
troscopy data in low fields B=l T-5.5 T together with B=0 T and B=7 T data for comparison, (b) 
Landau Levels in higher fields B=7 T to B=12 T clearly resolving LL with indices N=-4, -3, -2, -1, 
0. (c) Energy of the Landau Levels as function of sgn(N) yf\ N\ B. From the linear dependence, 
specific of Dirac fermions the Fermi velocity is extracted. 

We note that the LL at the Dirac point (DP) where N=0, is ~ 200 meV above the Fermi 
level, indicating that the system is hole doped even though the gate is at ground potential. 
The corresponding carrier density, n=2 x 10 12 cm -2 is found using: \Ep — Ej)\=%VF^[wri. 
This residual doping level, typical to graphene samples on SiC>2, can be attributed to charge 
impurities trapped between graphene and the SiC>2 or in the oxide. The energy of the DP 
relative to the Fermi level exhibits a small oscillation with field which, as discussed below, 
reflects the pinning of the Fermi level within each LL as the filling factor changes. This 
effect is much more pronounced when the filling factor is varied through gating as shown in 
Fig. i 

The absence of LLs in earlier STS measurements of graphene on SiC>2 despite evidence 
of their existence inferred by less direct probes^^ was puzzling^. To shed light on this 
question, we performed the same measurement on samples in which the SiC>2 dielectric 
was not chlorinated. The higher concentration of trapped charges in the unchlorinated 
substrates leads to a stronger random potential which is responsible for charge inhomogeneity 
and electron hole puddles in graphene. For these samples, even in those cases were peaks 
resembling LLs could be observed, their field and level index dependence did not follow the 
scaling expected for single layer graphene. Furthermore, similarly to earlier reports^, we 
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encountered gating hysteresis where sweeping the gate back and forth did not give identical 
spectra. We conclude that in order to access the intrinsic properties of the charge carriers in 
graphene it is essential to minimize the effect of trapped charges by using an appropriately 
treated dielectric substrate. 

The V shaped density of states in graphene gives rise to an ambipolar electric field effect 
which allows tuning the charge carrier density from electrons to holes by applying a gate 
voltage. Thus, for the samples studied here which are initially hole doped it is possible, by 
applying a positive gate, to bring the DP towards E F , as illustrated in Fig. md). In the 
STS experiment by changing the sample tip bias we probe the occupied and unoccupied 
states above and below E^. Therefore, when presenting the data it is convenient to keep 
the energy origin at Ep so that the energy of the DP varies with respect to it. To study the 
effect of gating on the LLs we record the differential conductance spectra, dI/dV(E), at a 
fixed value of the magnetic field for a sequence of gate voltages. Our results for B=12 T are 
summarized as a map of the spectra versus gate voltage in Fig. |4](a). Each vertical line in 
the map corresponds to the spectrum taken at a particular gate voltage. At Vg=-15 V, the 
DP (N=0) is at ~ E=240 meV, followed by a sequence of clearly resolved LL N = -1, -2, 
-3 etc. We vary the gate voltage between: -15 V < Vcate <+43 V, corresponding to carrier 
densities: 3xl0 12 cm" 2 > n > -0.5 xlO 12 cm -2 . 

Several features in this data stand out. Varying the carrier density through gating is 
accompanied by pinning of E^ to each LL as it is filled and followed by a jump of E^ to 
the next LL when a level is full. Qualitatively, one can understand the step-like features as 
follows: The LL spectrum consists of peaks where the DOS is large separated by minima 
with low DOS. It takes a large change in the charge carrier density to fill the higher DOS 
regions, resulting in plateaus where the Fermi level is "pinned" to a particular Landau level. 
In contrast, filling the low DOS region in between the LLs does not require a large change in 
carrier density therefore the jumps (changes in slope) in between plateaus. The broader the 
Landau levels are, the larger the slope in the plateau region and the less abrupt the jumps in 
between. A similar effect was previously observed in very high mobility GaAs samples using 
time domain capacitance spectroscopy (TDCS) 25 . Unlike the case of the 2DEG, due to the 
fact that in graphene the LL are not equally spaced, the largest jump from N=-l to N=0 
is followed by successively smaller jumps for higher index levels. Moreover, the broader the 
levels, the less abrupt the jumps, meaning that in disordered samples, this pinning effect is 
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FIG. 4: (a) Measured STS as function of gate voltage Vg (carrier density n). Each vertical line 
in the map is a spectrum at a particular gate. The bright stripes correspond to Landau Levels as 
indicated by N=0,± 1,± 2. (b) Simulation of the evolution of the LL spectrum with gate voltage, 
assuming Lorenzian line shapes of width ~ 30 meV, for all levels and v^=(1.16 ± 0.02) x 10 6 m/s. 
The latter was obtained from fitting the LLs at Vg=0 V , and C=100 aF/im -2 for the capacitance 
between the graphene and the Si back gate. This data and that in Fig. 3 were taken in different 
parts of the sample, therefore the small variation in Fermi velocity, (c) Gate voltage dependence of 
Vfi the average slope of the dispersion. The Dirac point, at Vc= 35 V, is marked, (d) Illustration 
of the effect of gating on the Fermi level. 
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smeared out. 

To analyze the data we model the Landau levels as Lorenzians with equal width w 30 
meV from which we numerically calculate the chemical potential as a function of carrier 
density. The result of this calculation shown in Fig. |4|b) is consistent with the experimental 
data. We note that as N=0 is brought closer to the Fermi level it becomes sharper and 
better defined as expected from theoretical considerations of minimal scattering at the Fermi 
energy^. 

For a linear dispersion (Dirac cone) one can obtain the Fermi velocity by plotting the LL 
energies against sgn(N) ■ \/\N\ B and fitting to Equation [IJ However when the dispersion 
is not strictly linear, as expected in the presence of electron-electron interactions, this pro- 
cedure gives the average slope of the dispersion,^. Using this procedure for a series of gate 
voltages we found a substantial increase in Vpupon approaching the Dirac point, Fig. (i[c), 
consistent with the expected interaction-induced renormalization of the Dirac con d 26 * 27 l 

By using STM/STS in magnetic field we have shown here that for graphene supported 
on chlorinated Si0 2 substrates, the substrate disorder is weak enough to allow observation 
of quantized LLs already at moderate fields. The ability to control the doping level in this 
system made it possible to observe the interaction-induced renormalization of the Dirac cone 
at low carrier densities. These results show that intrinsic phenomena reflecting correlation 
effects between the massless Dirac fermions are observable and become apparent even in 
non-suspended samples for sufficiently clean substrates. 

In a recent post Elias et alP^ report up to a three-fold interaction-induced increase in 
the Fermi velocity at low carrier densities obtained from the temperature dependence of the 
SdH in suspended graphene devices. 
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